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Abstract: 

A parameterized reaction mechanism has been 
created to examine the interactions between isoprene 
and other tropospheric gas-phase chemicals. Tests 

of the parameterization have shown that its results 
match those of a more complex reaction set to a high 
degree of accuracy. Comparisons between test runs 
have shown that the presence of isoprene at the 
start of a six day interval can enhance later ozone 
concentrations by as much as twenty-nine percent 
The test cases used no input fluxes beyond the 
initial time, implying that a single input of a 
biogenic hydrocarbon to an airmass can alter its 
ozone chemistry over a time scale on the order of a 
week. 

Introduction: 

Various studies (cf. Jacob and Wofsy, 1988, 
Trainer et al., 1987) have suggested that fluxes of 
terpenes such as isoprene may be responsible for 

enhanced ozone concentrations in the lower 

troposphere. The full reaction mechanisms for these 
gases are too complex for use in chemical transport 
models, due to the memory space required to store 
all species, and the computation time needed to 
predict future concentrations. The present work 
discusses the parameterization of one of these 
mechanisms (that of isoprene) into a simplified 

form. The effect of an initial input of isoprene on 
later ozone chemistry was studied as part of the 
parameterization process. 

Constructing the Parameterization. 

In the following discussion, reference will be 

made to three reaction sets; "no isoprene", 
"parameterization" and "full isoprene". The first 
is a simplified tropospheric reaction set, following 
Lurmann et al. (1986) and DeMore et al (1990). The 
"full isoprene" set includes the no isoprene 

reactions and the detailed isoprene system of 

Lurmann et al (1986), and incorporates some of the 
changes suggested by Jacob and Wofsy (1988), The 
"parameterization" is the simplified isoprene 
mechanism, designed to have the same effect on the 
simple troposphere as the full isoprene set, yet 
using fewer species to accomplish this goal. The 
parameterized reaction mechanism and its rates are 
given in Appendix I. 

The parameterization was constructed by 

eliminating unimportant reaction pathways and 

lumping species with similar chemistry together. In 
summary, the changes were: 

(1) Criegee - type radicals MVKO, MAOO, and MCRG 
were lumped together as "CREB". The products of the 
oxidation of CREB by NO were determined by adding 
the three oxidation reactions of CREB’s components. 

(2) Radicals MAN 2 and MVN 2 (produced via N0 3 
oxidation of methacrolein and methyl vinyl ketone, 
respectively) were lumped together as "MBN 2 ". The 


product coefficients for the subsequent reaction of 
MBN 2 with NO were derived by examining the amount of 
MBN 2 resulting from each of the methacrolein and 
methyl vinyl ketone paths. 

(3) Radicals MR0 2 and VR0 2 (produced via OH 
reactions with MACR and MVK, respectively) were 
lumped together as "BR0 2 ", with product species and 
reaction coefficients adjusted as in step (2). 

(4) The cycles for the natural acetyl nitrate 
species IPAN and MPAN were lumped together into the 
PAN cycle. The products leaving the cycle via 
oxidation of the lumped carbonate by NO are 
dependant on the concentrations of MACR and MVK. 
With this change, CH 3 CHO and HOCH 2 CHO are lumped as 
a single higher aldehyde, CH 3 C0 3 and MA0 3 are 
similarly lumped, and the total acetyl- nitrate is 
the sum of PAN, IPAN, and MPAN concentrations. The 
reaction rates used were those of the normal PAN 
cycle, but the products and their coefficients were 
changed to reflect the new sources of total PAN. 

(5) Product species IPN 4 , MGLO, pyruvic acid and 
methyl acrylic acid were dropped from the 
parameterization, as were the reactions for the 
production of ozonides. 

Numerical Tests of the Parameterization: 

Nine test simulations were performed, each 

simulation spanning six days of diumally varying 
photochemistry. No fluxes or deposition were 
allowed after the initial start-up of the box model. 
The initial conditions for these tests are given in 
Table 1, Most of the chemicals had the same initial 
concentration regardless of the test. The only 
species with different concentrations at the initial 
time were N0 2 , NO, 0 3 , PAN, and isoprene. NO x , 
ozone and PAN started at three levels, representing 
background, rural and polluted urban air. Isoprene 
was set at three levels. The combinations of the 
three NO^ scenarios with the three isoprene levels 
resulted in the nine test cases. Similarly, the 
level of NO x and isoprene are the only possible 

causes of the differences between the simulations. 

Reaction rates were calculated assuming sea-level 
pressures at the equator close to equinox. 

Photochemical rates were calculated using a detailed 
one-dimensional model (Henderson et al., 1987). The 
diurnal temperature profile and an example 

photochemical rate are given in Figure 1. 

Simulations began at noon on the first day, with 

the noon-to-noon diurnal profile in rates repeating 

during the subsequent five days. 

Two types of analysis were made with the 

resulting data. The effect of isoprene reactions on 
the model atmosphere were shown by comparing no 
isoprene and full isoprene simulations. The 
parameterization results were compared to the full 
isoprene results in order to gauge the usefulness of 
the former in simulating isoprene chemistry. 
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Figure 1. Temperature profile and J N02 . 



Only a few species will be shown and discussed, 
due to the large amount of data resulting from these 
simulations and the limited space for this 

presentation. Researchers interested in obtaining 
the data files or plots of the results should 

contact the authors at the above address. 

The 0 3 results for each test are given in Fig. 2. 
Each graph shows the concentrations (ppbv) as 
functions of time, and are labled according to the 
initial concentration scenario giving rise to the 
data. As can be seen from the graphs, the ozone 

concentration is enhanced by the addition of 

isoprene if the initial NO x level is high or medium, 
and is depleted if the initial NO x is low. The 
largest 0 3 increase is for the case (high initial 

NO x , high initial isoprene), in which the ozone 
maximum is 25 ppbv higher after the addition of 
isoprene (maximum of 83.8 ppbv without isoprene, and 


108.3 ppbv with initial isoprene), an enhancement of 
29 % . The parameterization’s 0 3 values were within 
a few percent of those from the full isoprene 
mechanism for all test runs. The parameterization 
would appear to be suitable for 0 3 simulations, with 
the caveat that ozone may be overpredicted by up to 
5 ppbv. 

The cause of the 0 3 enhancements are explained 
via the relative steady state between 0 3 and NO x . 
Steady state ozone concentrations are roughly 

proportional to the ratio of N0 2 to NO (0 3 a 
Ji«NO ? /(k 14 «NO)). Reactions that bias this NO x 
ratio in favour of N0 2 therefore result in higher 
ozone concentrations. 

Two processes cause such a bias. The first is 
the reaction of NO with R0 2 radicals to produce N0 2 . 
These radicals (RI0 2 , IN0 2 , VR0 2 , MVN 2 , MR0 2 , and 
MAN 2 ) originate from isoprene’s oxidation by OH and 
N0 3 . The reaction rates for the first hour of the 
"high NO x - high Isoprene” scenario have shown that 
43% of the net (positive) rate of change of N0 2 is 
due to production of N0 2 via the R0 2 reactions. The 
addition of isoprene causes the N0 2 /N0 ratio to move 
from 8.13 to 9.09. 

This process is short-lived due to the rapid 
oxidation of isoprene and its product ketones within 
the first few hours of the simulation (example 
isoprene concentrations are given in Figure 3). 
This removes the source term for the R0 2 radicals. 
The second increase in 0 3 , between hours 18 and 24 
(6 am and noon) is due to the release of N0 2 from 
the cycles of peroxyacetylnitrate (PAN) , MPAN, and 








IPAN. These species have very similar reaction 
chemistry (hence their combination in the 
parameterization); production is via reaction of a 
peroxyacetyl radical with N0 2 , and loss is via 
thermal dissociation of the compound back into N0 2 
and the radical. The colder temperatures during the 
night (hours 6 to 18) thus resulted in a buildup of 
the PAN-type species. The increase of temperatures 
on the dawn of the second day caused the release of 
N0 2 stored during the night as PAN, MPAN and IPAN. 
These processes have a magnitude sufficient to 
account for the entire net rate of change of N0 2 at 
11:30 am on the second day. 

The enhancement of 0 3 due to PAN-type compounds 
from isoprene oxidation has important implications 
with regards to ozone production far from the source 
region. For example, NO x produced in a region of 

biomass burning (a high NO x , high isoprene source 
region), if carried aloft a few kilometres (a 
temperature contrast equivalent to the day/night one 
used here) could be stored as PAN and related 
organic compounds. Later heating of the airmass far 
downwind from the source would result in strongly 
enhanced ozone. Such a mechanism may be responsible 
for the "Ozone high" in the troposphere over the 
South Atlantic, mentioned in other papers in these 
Proceedings. 



isoprene mechanism. 

The mechanism responsible for ozone depletion in 
the "low NO " simulations is that of removal by 
isoprene. Initial NO x levels are sufficiently low 

so that the (R0 2 + NO) ozone source fails to 
compensate for the loss of ozone due to isoprene 
oxidation. During the first nightfall, isoprene 

removal of ozone dominates the other loss terms by 
about two orders of magnitude. Thereafter the curve 
describing ozone concentrations closely follows that 
of isoprene, with a gradual decrease in both 

species. 

Total lumped PAN’s for high NO x - high isoprene 
are given in Figure 4. Total PAN has increased by 
about an order of magnitude for the high and medium 

NOx scenarios. The increase is smaller for the low 

NOx case, and almost non-existant for the case of 

low-NOx, high-isoprene. The enhancement of total 
PAN’s is due to the addition of the biogenic MPAN 
and IPAN to peroxyacetylnitrate. These species 

increase in concentration until the middle of the 
second day, when the concentrations of their ketone 
precursors have become depleted. The enhanced PAN 
cycle is one of the primary causes of the ozone 
enhancement noted above. 


Methane, ethane and propane all have a 
"staircase" pattern of decreasing concentration 
(ethane is given as an example, Figure 5). The most 
rapid removal (steep parts of the steps) occurs 
during the day, when the oxidant, OH, is highest in 
concentration. The differences in the rate of 
alkane removal between isoprene and no isoprene 
cases thus reflect differences in OH concentrations. 
Ethane removal is enhanced by the addition of 
isoprene in the high and medium NO x scenarios (ie. 
higher OH). In the low NOx scenario, the addition 
of isoprene results in a decrease in the ethane 
removal rate (ie. lower OH). The OH values are a 
direct reflection of the NO x changes discussed 
above. In the high NO x - high isoprene case, 
daytime OH concentrations are enhanced by higher 0 3 
and thus higher 0( ] D) values. In the low NO x - high 
isoprene case, daytime OH is inhibited by removal by 
the as yet undepleted isoprene. At hour 24 for this 
scenario, the loss rate of OH due to isoprene alone 
is 4.2 times the combined loss rates due to CO and 



Figure 4. Total PAN PAN + MPAN + [PAN) 
Concentrations. 



CH 4 oxidation. When combined with the effects of 
its ketone products, methacrolein and methyl vinyl 
ketone, the isoprene-caused OH loss rate is 5.2 
times the combined CO and CH 4 loss rate. In regions 
where the isoprene concentration is high (eg. close 
to the source) it will be the main chemical 
determining OH concentrations. 

Formic acid values are given in Figure 5 (again, 
high NO x - high isoprene is used as an example). 
The production of formic acid is greatly enhanced by 
the introduction of isoprene for all test cases. 
'Phis is due to the addition of the creigee biradical 
CRF.A (CH 2 0 2 ) as a source of formaldehyde. 






Conclusions: 

The isoprene parameterization derived in this 
work has been compared to a more detailed reaction 
mechanism, and has been found suitable for use in 
applications in which the number of chemicals are 
limited by computational memory and processing time. 
The parameterization uses 14 less species than the 
full mechanism, yet produces essentially the same 
results for the most important tropospheric gases. 

The effect of isoprene on tropospheric chemistry 
depends on both the initial NO x and isoprene levels. 
Ozone was found to either increase or decrease in 
concentration depending on the initial NO x values. 
Increases were due to conversion of NO to N0 2 , with 
the causes of this conversion being R0 2 reactions 
with NO x and release of N0 2 from PAN -type 
compounds. Similar effects were shown for the other 
species studied. The lack of deposition or sources 
after the initial time showed that the injection of 
isoprene into an air mass can have a long term 
effect on the chemistry of that air mass, long after 
the precursor hydrocarbon has been totally oxidized. 

Table 1: Initial Concentrations 


Sped** 


Concentration (ppbv) 


Concentration (molecule* errr 1 ) 

NO, 


0.1 / 3.0 / 8.0 


0.25 / 7.4 / 20.0 (xlO 10 ) 

NO 


0.02 / 1.0 / 1.8 


0.05 / 2.5 / 4.4 (xlO 10 ) 



20. / 50. / 80. 


50. / 125. / 200. (xlO 10 ) 

NO, 


0,04 pptv 


10 s 

N,Oj 


0.04 pptv 


10* 

PAN 


0.01 / 0.3 / 8.0 


0.25 / 7 5 / 20. (xlO 9 ) 

HNO, 


02 


5.0 (xlO 1 *) 

HNO, 

20 


5.0 (xlO 10 ) 

hno 4 


1.0 pptv 


2.5(xl0 2 ) 

CH, 


1.67 ppmv 


4.185 (xlO 11 ) 

HCHO 


2.0 


5.0 (xlO 10 ) 

MCHO 

0.5 


1.25 (xl0 lc> ) 

MCO, 

; 004 pptv 


10* 

MOOH 


02 


5.0 (xlO*) 

MO, 


0.04 pprv 


10* 

MO 


0.04 pptv 


10* 

CO 


200 


5.0 (xlO 12 ) 

CO, 


343 ppmv 


8 575 (xlO 15 ) 

CjH* 


0.8 


I 2.0 (xlO 10 ) 

C,H, 


0.2 


5.0 (xl0 lflr ) 

ETO, 


0.04 pptv 


10* 

H,0, 


0.5 


1.25 (xl0’°) 

HO, 


0.04 pptv ' 

10* 

OH 


0.04 pptv 


to* 

FORM 


0.1 


2.5 (xlO 10 ) 

AHO, 


0.04 pptv 


10* 

0 ! D 


4.0 (xlO 10 ) 


10 1 

LOSS 


4.0 (* lO s0 ) 


10' 

HjO 


15 % 


3 829 (xl0 |T ) 

Oj 1 



5 131 (xlO 1 *) 

M 



2.45 (xl0 l *l 


Full boprwe Miriwlw InfeM Cmmkrnt 


Spade* 

" 

Concentration (ppbv) 

— 

Concentration (molecules cm* 2 ) 

ISOP 

MACR 

MVK 

RIO, 

INO, 

MGGY 

IPAN 

MPAN 

HAC 

PYRU 


03 / 3.1 / 8.2 
0.2 
0.2 

0.04 pptv 
0.04 pptv 
4 pptv 
4 pptv 
4 pptv 
0.01 

0.04 ppcv 
(all remaining 
species, same 
as PYRU) 


1.23 / 7.30 / 20.0 (XlO 40 ) 

5.0 (xlO 1 *) 

5.0 (xlO 10 ) 

10 * 

10 * 

10 * 

10 * 

10 * 

2.5 (xlO*) 

10 * 

(all remaining species; 
MAAC VRO,, MAN,, MVN, t 
MROj, IPN 4 , CREA, MVKO, 
MAOO, MCRO, MGLO. HACO. 

MAO,, and OTTO, start 

at 10* nwlsmlss car*). 


ISOP, MAC*, MVK. RKX INO. MOO Y . CUA an seme m above. 
The redial* CUB. MBN» and 8*0, all start at 10* molecules 
car* The lumped PAN ha* initial concemndoae of 43x10*, 
7 .7x10*, end 2 .02x10* molecule* cm’, far low, medium end high 
NCL scenario*. The lumped MCHO concentration is 1.27SxlO>* 
molecules car* 


Appendix I : Reaction Mechanisms 


TYPE SOURCE 


1 NO, 4 hv 

-* NO 

4 0, 

1 

a 


2 0 , 4 hv 

4 0( S D) 

♦ °I 

1 

a 


3 HA + hv 

4 2 OH 


1 

a 


4 HCHO + hv 

■* H, 

4 CO 

1 

a 


5 HCHO + hv 

4 2 HO, 

4 CO 

1 

a 


6 MOOH + hv 

4 MO 

4 OH 

1 

a 


7 N,0, + hv 

4 NO, 

4 NO, 

1 

a 


8 HNO, + hv 

4 OH 

4 NO, 

1 

a 


9 MCHO + hv 

4 MO, 

4 HO, 4 CO l 

a 


10 NO, + hv 

4 NO 

4 O, 

1 

a 


11 NO, 4 hv 

4 NO, 

4 0, 

1 

a 


12 OOD) 4 M 

40, 

4 M 

special function 4 

a 


13 O0D> 4 H z O 

4 2 OH 


2.2E-10 3 

a 


14 O, 4 NO 

4 NO, 

40, 

2.0E-12 1.4E403 2 

a 


15 O, * OH 

4 HO, 

40, 

L.6E-12 9.4E402 2 

a 


16 O, 4 HO, 

4 OH 

4 2 0, 

LIE-14 5.0E402 2 

a 


17 O, + NO, 

4 NO, 

4 0 , 

1.2E-13 2.45E403 2 

a 


18 OH + H,0, 

4 HO, 

4 H,0 

2.9E-12 1.6E402 2 

a 


19 OH ♦ NO 

4 HNO, 


special function 4 

b 


20 OH + NO, 

4 HNO, 


special funedoo 4 

a 


21 OH + HNO, 

4 NO, 

4 H,0 

special funedoo 4 

a 


22 OH 4 CO 

4 HO, 

4 CO, 

special funedoo 4 

a 


23 OH + CH 4 

4 MO, 

4 H,0 

2.3E-12 1.7E+03 2 

a 


24 OH 4 HCHO 4 HO, 

4 CO 4 H-0 l.OE-ll 3 

a 


25 OH + MOOH MO, 

4 H,0 

3.8E-12 *2.0E+O2 2 

a 


26 OH + MCHO 4 MCO, 

4 H,0 

6.9E-12 -2.5E+02 2 

b 


27 HO, 4 NO 

4 NO, 

4 OH 

3.7E-12 -2.4E+02 2 

a 


28 HO, + OH 

4 H,0 

40, 

4.8E-1 1 -2.5E+02 2 

• 


29 2 HO, 

4 H,0, 

+ O, 

2.3E-13 -6.0E+02 2 

a 


30 HO, + MO, 

4 MOOH 

40, 

3.3E-13 -8.0E+02 2 

a 


31 HO, + MCO, 

4 LOSS 


3.0E-12 3 

b 


32 HO, + NO, 

4 hno 4 


special funedoo 4 

b 


33 NO, + NO, 

4 N,0, 


special function 4 

a 


34 NO, + NO, 

4 NO 

4 NO, 4 O, 

2.5E-14 1.23E403 2 

b 


35 N.O, 

4 NO, 

4 NO, 

special function 4 

b 


36 NO, 4 NO 4 

2 NO, 


8.0E-12 -2.5E402 

2 

b 

37 NO, + HCHO 4 

HNO, 

4 HO, 4 

CO 3.2E-16 

3 

b 

38 NO 4 MO, •* 

NO, 

4 MO 

4.2E-12 -1.8E+02 

2 

b 

39 NO, > MCO, 4 

PAN 


4.7E-12 

3 

b 

40 NO + MCO, 4 

MO, 

- CO, 4 

NO, 4.2E-12 -I. 8 E 4 O 2 

2 

b 

41 PAN 

MCO, 

♦ NO, 

1.95E416 9.0E-402 

2 

a 

42 MO >0, 

HCHO 

4 HO, 

3. 94 E- 14 9.0E402 

2 

a 

43 N,0, 4 H,0 4 

2 HNO, 


1.3E-21 

3 

b 

44 CHi OH 

ETO, 

+ H,0 

1.7E-1 1 1 232E403 

2 

b 

45 ETO, *■ NO 4 

MCHO 

4 HO, 4 

NO, 4.2E-12 -L8E-402 

2 

b 

46 2 ETO, -» 

16 MCHO 

4 1.2 HO, 

5.0E-14 

3 

b 

47 ETO, + HO, 4 

LOSS 


3 OE-12 

3 

b 

48 C,h/ + OH 

ETO, 


1 18E-11 6.79E402 

2 

b 

49 MCHO ♦ NO, - 

MCO, 

4 HNO, 

1.4E-12 1.9E403 

2 

a 

50 HNO, 4 

NO, 

4 HO, 

special function 

4 

b 

51 HNO, + hv 

NO 

4 OH 

0 205 x /(NO,) 

1 

b 

52 HCHO + HO, 4 

AHO, 


1.0E-14 

3 

b 

53 AHO, + NO 4 

FORM 

4 HO, 4 

NO, 4 2E-12 -1.8E402 

2 

b 

54 AHO, + HO, 4 

FORM 

4 H,0 4 

O, 2 OH-12 

3 

b 

55 2 AHO, 

2 FORM 

4 2 HO, 4 

2 O, 1.0E-13 

3 

b 

56 FORM 4 OH 4 

HO, 

4 H,0 4 

CO, 3.2E-13 

3 

b 


Comments: 

( 1 } Photolysis rate* from detailed model result* at 0 km (see text). 

(2) Temperature de pendant rates, k ■ k, x exp4 -kjT), T -semp. in K. 


Q.-A 


(3) Constant rates. 






+ 0.33 MACR 4 0.33 MOGY 


(4) Special functions for given reaction*: 

(12) The reaction O0D) 4 M 4 CX 4 M (the product* an actually 
OOP) M, but the conversion to O, la assumed to be 
"instantaneous") is a combination of the reactions 

CK l D) ♦ N, 4 OOP) + N a . k-l.8E-tl • e ll0fr , and 

O0D) 4 Q, 4 OOP) + O, . k-3.2E-Ll • e™ 

If 80 % of M is Nj, and 20 % is Oj, then the net reaction is 
O0D) 4 M 4 O, + M. and 


(19) k 


- 1.44E-1 l«e 
A*T a * 


HOT 


1 4 A*' 

1.0E-22. B 




C»T U 
3.3, C ■ 


[M] 


4 6.4E-I2«c 

0.6 


,*YT 


(20) k ■ 


A* ( M ] 


1 4 A* 




. 9.0E-09, D 

fl + (lo«,o 
► 0.6 ^ 


( lt ( l 0!ir A^Mi)y 
UiL ) 'ji 


1.0 

A 


A - 2.6E-30 • (T/300)' . B - 2.4E-U • (T/300) 


(21) k 


a C *[M] 

A 4 - — , A 


1 4 C* [M] 

B 

4. IE- 16 • exp(1440/D, C 


- 7.2E-15 • exp ( 783 / T ), 
1.9E-33 • exp (723 / T ). 


(22) k - 1.50E-13 • (1 + 0.6»[M]*T* 1.36E-16 / 1.013E+06) 
f 6.9E-33*rM]*e 1 007,T ) 

(32) k ll . 4.86E- 1 2*(M] U M J 

(33) (same formula as reaction 20), 

A - 2.2E-30 * (T/300)'*" 

(35) k - k„* 7.5E+26 *(300/T) OJ 
(50) k - ( 4.9E-06«(M]*exp(-l0015/T) ) / (1 

In the above, [M] is the total number density in molecules cm 3 , and T 
is the temperature in degrees Kelvin. 

Sources: (a) DeMore et al, 1990, (b) Lurmann et aL, 1986 


I 1.5E-12 • a/300r 

-UMOT 


- 4.86E-12*[M] ) 


Parameterization Reactions: 


57 ISOP 

+ 

O, 

4 

0.5 HCHO 

4 0.2 MVK 

4 

0.30 MACR 






4 

0.2 CREA 

+ 0.06 HO, 

4 

0.5 CREB 

7.0E-15 1.9E403 

2 

58 MVK 

4 

O, 

4 

0.5 HCHO 

4 0.2 CREA 

4 

0.21 HO, 






4 

0.2 CREB 

+ 0.15 MCHO 

4 

0.5 MGGY 






4 

0.15 MCO, 




4.0E-15 2 0E403 

2 

59 MACR 

4 

O, 

4 

0.5 HCHO 

4 0.2 CREA 

4 

0.21 HO, 






4 

0.15 MO! 

4 0.5 MGGY 



4.4E-15 2.5E403 

2 

60 ISOP 

♦ 

OH 

4 

RIOj 




1.5E-11 -5.0E+O2 

2 

61 RIO! 

+ 

NO 

4 

0.9 NO, 

4 0.45 MVK 

4 

0.45 MACR 






4 

0.9 HO, 

4 0.9 HCHO 



4 2E-12 -1.8E402 2 

62 ENOj 

4 

NO 

4 

2.0 NO, 

4 HCHO 

4 

0.5 MVK 






4 

0.5 MACR 




42E-12 -1 8E+02 

2 

63 tNQ, 

4 

NO, 

4 

LOSS 




4.2E-13 -1.8E+02 

2 

64 [SOP 

4 

NO, 

4 

INOj 




3.23E-13 

3 

65 CREA 

4 

H,0 

4 

FORM 

4 H,0 



4.0E-18 

3 

66 CREA 

4 

NO 

4 

HCHO 

4 NO, 



7.0E-12 

3 

67 CREA 

4 

no 2 

4 

HCHO 

4 NO, 



7.0E-13 

3 

68 CREB 

4 

HjO 

4 

LOSS 

4 H,0 



4.0E-18 

3 

69 CREB 

4 

NO 

4 

0 33 MVK 

4 0.33 MACR 

4 0.33 MGGY 






4 

NO, 




4.2E-12 -1.8E+02 

2 


70 CREB ♦ NO, 4 0.33 MVK 



4 NO, 



4.2E-13 -1.8E402 2 

71 RIO, 

4 HO, LOSS 



3.0E-12 3 

72 INO, 

4 HO, 4 LOSS 



3.0E-12 3 

73 MACR 

4 NO, 4 MBN, 



6.7E-15 3 

74 MVK 

4 NO, 4 MBN, 



6.0E-14 3 

75 MBN, 

4 NO 4 2.0 NO, 

4 HCHO 

4 0.45 HO, 



4 0.45 MCO, 

4 0.55 MGGY 


4.2E-12 -1.8E402 2 

76 MACR 

4 NO, 4 MCO, 

4 HNO, 


3.3E-15 3 

77 MGGY 

4 MCO, 

4 HO, 


0.15 xJ(NO,) 

78 MGGY 4 OH 4 MCO, 



1.7E-11 3 

79 MACR 

4 OH 4 BRO, 



3.86E-12 -5.OE402 2 

80 MACR 

4 OH 4 MCO, 



1.02E-1 1 

81 MVK 

4 OH 4 BRO, 



3.0E-12 -5.0E4O2 2 

82 BRO 2 

4 HO, 4 LOSS 



3.0E-12 3 

83 BRO, 

4 NO 4 0.9 NO, 

4 0.625 HO, 

4 0.625 HCHO 



4 0.623 MGGY 

4 0.275 MCO, 

4 0.275 MCHO 4.2E-12 -1.8E402 2 

84 MBN, 

4 HO, 4 LOSS 



3.0E-12 3 


Comments: 

( 1 ) 


The first 56 reactions 


those of the "No Isoprene" set. 


with the exccodon of reactioo 40, which hss been aliened' to reflea 
1 PAN cycle. Reaction 40 is replaced with: 


the lumped 

MCO, + NO 


t MO, + s COj + (a+3b) NO, 4 b HO, 4 b MGOY - 2b NO. 


k , 4 (MACRXNO, ) ♦ k, ,(MACR(OH) 

'(k^tNOjl+k,, (OH))(MACR) 4 O 1 * t ( N0 » > +k i ♦ < >) C MCHO) ♦ 0 15 k , , <MVK)(0,)j 

4(0 45 k M (MBN,) 4 0 2 7 3 k,, ( BRO, )) (NO) 4 (k T , 4 k t| (OH)) (MGOY > 

and s - 1,0 - b. In the parameterization, "MCO,", "MCHO". and "PAN" 
refer to the lumped peroxyacetyl radical, lumped higher aldehyde and 
lumped acetylnitate species, respectively. 

(2) Temperature dependant rates; k-k 1 *exp(-k v >T), T - 
temperature in degrees K. 

(3) Constant rates. 
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